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We study the 2-flavour lattice Schwinger model: QED in D = 2 with two fermion species of identical mass. In 
the simulation we are using Wilson fermions where chiral symmetry is explicitly broken. Since there is no known 
simple order parameter it is non-trivial to identify the critical line of the chiral phase transition. We therefore 
need to find observables which allow an identification of a possible restoration of chiral symmetry. We utilize the 
PCAC-relations in order to identify the critical coupling, where chiral symmetry is restored. 



1. INTRODUCTION 

The 2-flavour Schwinger model is a 2D theory 
of two fermion species and photons. In the mass- 
less continuum model the SU{2)a flavour symme- 
try is broken by a mechanism mimicking sponta- 
neous symmetry breaking. Due to the Mermin- 
Wagner-Coleman theorem ^ in this situation 
there is no spontaneous symmetry breaking and 
thus there are no Goldstone bosons, nevertheless 
massless pions are expected from continuum stud- 
ies (for a recent study cf. |^). 

The naive lattice formulation for fermions 
shows the notorious doubling problem. In Wil- 
son's action a term is added giving the doublers 
a mass 0{\/a) decoupling them in the contin- 
uum limit. Independent of the possible sponta- 
neous breaking of chiral symmetry (which does 
not occur in 2D but is expected in e.g. QCD4) 
we expect that the chiral symmetry of the ac- 
tion itself is restored at a critical fermion coupling 
Kc{(3). However, since the Wilson term breaks 
chiral symmetry explicitly we have no known sim- 
ple order parameter, which allows one to iden- 
tify that position. Following earlier proposals [|j 
it has recently been suggested to utilize PCAC- 
relations with Schrodinger functional methods for 
that purpose Q ; this turned out to be quite pow- 
erful in the context of QCD4 Here we report 
on results for QED2 with 2 species of fermions, 
where we also employed PCAC-relations in order 
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to determine the critical fermion coupling line. 

The lattice action for the massive 2-flavour 
Schwinger model with Wilson fermions is 

S[U,ij,^P]^Sg[U] + Sf[U,ij,^], (1) 

where Sg is the compact plaquette gauge action 
(coupling (3), ijj = (m, d) is a flavour-doublet and 



Sf[U,Tp,1p] = 2^ VnV'r! 

n 

-K ^ [■(/!„ (r - cr^) Un^f, Ipn+i, + h.c] 



(2) 



n,fj, 



The ipip term and the terms proportional to r 
(we choose r — 1) break chiral symmetry explic- 
itly. Strictly speaking the action has only the 
SU{2)v X U{l)v global symmetries. The main 
idea using the PCAC-concept is to find the points 
in the coupling space, where the explicit break- 
ing of the chiral symmetries SU{2)a x U{1)a 
vanishes. The fermions were simulated with the 
HMC-method; more details can be found in Q. 

2. WARD IDENTITIES 

Global symmetries of an action on the classi- 
cal level lead to conserved Noether-currents. In a 
QFT the symmetries may be spontaneously bro- 
ken or broken due to anomalies resulting from the 
functional integration. So the relations express- 
ing current conservation may change depending 
on the mechanism of symmetry violation. To de- 
rive these relations one considers global infinitesi- 
mal symmetry transformations of a path integral 

(0|O[f/,V5,V']|0) = I f -^rr,J.,j.,nnr J. J,.-S 



2 



0.20 



0.15 



"left 0.10 



0.05 



0.00 





0.35 



Figure 1. Mass plateaus at /3 = 2.0 for the 16 x 16 
lattice (10000 configurations). 

where O denotes an arbitrary operator. If the 
integration measure is invariant this leads to re- 
lations like 



(0|^O|0} - (0|O^S'|0} = 



(4) 



called Ward identities. Thus, the quantum analog 
to the classical conservation law has similar (or 
equal) form but is an operator identity. 

For the massive continuum Schwinger model 
with two flavours the corresponding identities for 
the global SU{2)v x C/(l)v symmetries which are 
exactly realized in the QFT are 

. SU{2)v d^r^ = 0, j^- = ^jll^^p,(5) 
. Uil)v d,f -0, 4,1^^1^ 

The SU{2)a x U{1)a symmetries are broken for 
massive fermions. The SU{2)a Ward identity 
(PCAC-relation) reads 



Figure 2. Values (and linear fits) obtained for 
TOeff for different k at a given [3 (from left to right: 
0.1, 1, 2, 4) on an 8x8 lattice (1000 configurations 
each). 

3. RESULTS 

The underlying concept is to utilize a lattice 
version of (0) in order to identify the value of 
Kc(/3) where SU{2)a holds. That point defines 
a theory with massless lattice fermions. The 
naively discretized continuum action 

S ^ Sg[U]+ TOefT ^ V'nV'n 



2 ^ 



SU{2)a 9^/°^ = 2mo7r",7r" 



leads to an SU{2)a lattice operator identity anal- 
ogous to (|^). In order to satisfy that relation 
non-trivially we choose to take expectation val- 
ues with regard to a pion source and obtain 

yn^y: 2 TOeff (0|<+ < |0) = 



J 



(7) 



The U{1)a symmetry is also broken by the 
anomaly 

• U {1)a ^2moi) ^^"^^ ® ^ ^ _f- anomaly , 



(8) 



E [(0| < - < |0)] (10) 



with 



i>n- 



tfj. 



(11) 
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Figure 3. Phase diagram for the 2-flavour model 
for the 16 x 16 lattice (dashed lines to guide the 
eye); the value for the 1-flavour model is from 0. 



The operators in the above relation are projected 
onto states with zero momentum by summing 
over time-slices. This relation enables us to mea- 
sure the mass parameter mgfF- Here we are only 
interested in the value Kc(/3), where Weff vanishes. 

The identity ( |lO| ) should hold independent of 
the time-separation between source and sink. 
In Q such a behaviour in the context of the 
Schrodinger functional is used to determine the 
value of improvement coefficients. In our case we 
find plateaus already as soon as the distance ex- 
tends beyond the overlap region of the operators 
(fig. 1). For the analysis we excluded the points 
up to distance 3 (for 8x8 lattice: 2). 

The remaining points represent a mass plateau 
giving the value of rh^fi and its error is obtained 
from the weighted average. The errors of the cor- 
relation functions were computed using the jack- 
knife method. 

For given (3 we determined the values of Weff 
for different values of k. According to the gen- 
eral idea chiral symmetry is restored where TOgflF 
vanishes. Interpolating the values of m^g we get 
the desired value of ndfi) (fig- 2). In the domain 
studied a linear fit is appropriate. The resulting 
phase diagram is given in fig. 3. 

Finally, fig. ^ indicates possible problems. 
When measuring TO^ff we simultaneously kept 
track of the topological charge i/ of the configura- 
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Figure 4. Dependence of Wefj (circles) on the 
topological charge i' {/3 — 6, lattice size 8x8, 
1000 configurations each; Uc denotes the number 
of configurations with non-zero v) . 



tions. In particular at larger /3 tunneling between 
different sectors becomes less frequent Our 
results demonstrate a distinct dependence of the 
measured m^g on v which may lead to problems 
in situations where a simulation is non-ergodic. 
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